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t"*^ , Abstract. After a brief review of the phenomena expected in cold dense quark mat- 

^S| ' ter, color superconductivity and color-flavor locking, we sketch some implications of 

f^ I recent developments in our understanding of cold dense quark matter for the physics 

Q\ . of compact stars. We give a more detailed summary of our recent work on crystalline 

^^ ' color superconductivity and the consequent realization that (some) pulsar glitches may 

originate in quark matter. 
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,__! . 1 Color Superconductivity and Color-Flavor Locking 

>: 

r^ ' Because QCD is asymptotically free, its high temperature and high baryon den- 

2-} ', sity phases are more simply and more appropriately described in terms of quarks 

^j. • and gluons as degrees of freedom, rather than hadrons. The chiral symmetry 

f^*) I breaking condensate which characterizes the vacuum melts away. At high tem- 

f^ ■ peratures, in the resulting quark-giuon plasma phase all of the symmetries of 

^P I the QCD Lagrangian are unbroken and the excitations have the quantum num- 

(~| ■ bers of quarks and gluons. At high densities, on the other hand, quarks form 

CIh! Cooper pairs and new condensates develop. The formation of such superconduct- 

jis ■ ing phases ||l|,||||,§,||,|| requires only weak attractive interactions; these phases 

qj I may nevertheless break chiral symmetry [g| and have excitations with the same 

r^ • quantum numbers as those in a confined phase [§|J3,@,y . These cold dense quark 

^ I matter phases may arise in the cores of neutron stars; understanding this region 
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of the QCD phase diagram requires an interplay between QCD and neutron star 
phenomenology. 

The relevant degrees of freedom in cold dense quark matter are those which 
involve quarks with momenta near the Fermi surface. At high density, where the 
quark number chemical potential /i (and hence the quark Fermi momentum) is 
large, the QCD gauge coupling g(/i) is small. However, because of the infinite de- 
generacy among pairs of quarks with equal and opposite momenta at the Fermi 
surface, even an arbitrarily weak attraction between quarks renders the Fermi 
surface unstable to the formation of a condensate of quark Cooper pairs. Creat- 
ing a pair costs no free energy at the Fermi surface and the attractive interaction 
results in a free energy benefit. Pairs of quarks cannot be color singlets, and in 
QCD with two flavors of massless quarks, they form in the (attractive) color 3 
channel in which the quarks in a Cooper pair are color-antisymmetric |H,pp,H| . 
The resulting condensate creates a gap A at the Fermi surfaces of quarks with 
two out of three colors, but quarks of the third color remain gapless. Five glu- 
ons get a Meissner mass by the Anderson-Higgs mechanism Q; a SU {2)coIot 
subgroup remains unbroken. The Cooper pairs are flavor singlets and no fla- 
vor symmetries are broken. There is also an unbroken global symmetry which 
plays the role of U{1)b- Thus, no global symmetries are broken in this 2SC 
phase. There must therefore be a phase transition between the hadronic and 
2SC phases at which chiral symmetry is restored. This phase transition is first 
order |^ ll,lj,|l3| since it involves a competition between chiral condensation 



and diquark condensation |Tl| , [l3|| . 

In QCD with three flavors of massless quarks, the Cooper pairs cannot be 
flavor singlets, and both color and flavor symmetries are necessarily broken. 
The symmetries of the phase which results have been analyzed in PJ^. The 
attractive channel favored by one-gluon exchange exhibits "color-flavor locking." 
A condensate of the form 

{i'Ti'f) a Z\e"'^^e"''^ (1) 

involving left-handed quarks alone, with a, (3 color indices and a, h flavor indices, 
locks SU{3)l flavor rotations to SU{3)co\oi'- the condensate is not symmetric 
under either alone, but is symmetric under the simultaneous SU{3)l+coIot: ro- 
tations. R A condensate involving right-handed quarks alone locks SU{3)ii flavor 
rotations to SU{3)co\oi- Because color is vectorial, the combined effect of the LL 
and RR condensates is to lock SU{3)l to SU{3)r, breaking chiral symmetry.n 
Thus, in quark matter with three massless quarks, the SU(3)coioi x SU{3)l x 



It turns out Iq] that condensation in the color 3 channel induces a condensate in 
the color 6 channel because this breaks no further symmetries M. The resulting 
condensates can be written in terms of ki and K2 where {iI'l'^'4'i^) ~ ki(5""5'^'' -|- 
K.25°'''S^°'. Here, the Kronecker S's lock color and flavor rotations. The pure color 3 
condensate (hi) has K2 ~ — ki. 

Once chiral symmetry is broken by color-flavor locking, there is no symmetry argu- 
ment precluding the existence of an ordinary chiral condensate. Indeed, instanton 
effects do induce a nonzero (qq) M, but this is a small effect ]l4|]. 
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SU{3)b.^U{1)b symmetry is broken down to the global diagonal SU{3)coior+L+R 
group. A gauged U{1) subgroup of the original symmetry group — a linear 
combination of one color generator and electromagnetism, which lives within 
SU{3)l X SU{3)r — also remains unbroken. All nine quarks have a gap. All eight 
gluons get a mass pjlq]. There are nine massless Nambu-Goldstone bosons. All 
the quarks, all the massive vector bosons, and all the Nambu-Goldstone bosons 
have integer charges under the unbroken gauged U{1) symmetry, which there- 
fore plays the role of electromagnetism. The CFL phase therefore has the same 
symmetries (and many similar non-universal features) as baryonic matter with a 
condensate of Cooper pairs of baryons . This raises the possibility that quark 
matter and baryonic matter may be continuously connected ||7|] . 

Nature chooses two light quarks and one middle-weight strange quark, rather 
than three degenerate quarks. A nonzero rus weakens condensates which involve 
pairing between light and strange quarks. The CFL phase requires (us) and (ds) 
condensates; such condensates which pair quarks with differing Fermi momenta 
can only exist if the resulting gaps are larger than of order ms/2/z, the difference 
between u and s Fermi momenta in the absence of pairing. This means that upon 
increasing nis at fixed n, one must find a first-order unlocking transition H,^: 
for larger nis only u and d quarks pair and the 2SC phase is obtained. For any 
m-s ^ oo, the CFL phase is the ground state at arbitrarily high density Q. 
For large values of rris, there is a 2SC interlude: as a function of increasing /i, 
one finds a first order phase transition at which hadronic matter is replaced by 
quark matter in the 2SC phase and a subsequent first order phase transition at 
a higher fi above which CFL quark matter takes over. For smaller values of tti^, 
the possibility of quark-hadron continuity [[7| arisesjj 

Much effort has gone into estimating the magnitude of the gaps in the 2SC 
and CFL phases, and the consequent critical temperature above which quark 
matter ceases to be superconducting. It would be ideal if this task were within 
the scope of lattice gauge theory. Unfortunately, lattice methods relying on 
importance sampling have to this point been rendered exponentially imprac- 
tical at nonzero baryon density by the complex action at nonzero /i. There are 
more sophisticated algorithms which have allowed the simulation of theories 
which are simpler than QCD but which have as severe a fermion sign prob- 
lem as that in QCD at nonzero chemical potential |l6[. This bodes well for 
the futuren Given the current absence of suitable lattice methods, the mag- 



Note that even if the strange and light quarks are not degenerate, the CFL phase 
may be continuous with a baryonic phase which is dense enough that the Fermi 
momenta of all the nucleons and hyperons are comparable; there must, however, be 
phase transition(s) between this hypernuclear phase and ordinary nuclear matter M. 
Note that quark pairing can be studied on the lattice in some models with four- 
fermion interactions and in two-color QCD |17|. The Nc — 2 case has also been 
studied analytically in Refs. l4,[18|; pairing in this theory is simpler to analyze because 
quark Cooper pairs are color singlets. The Nc —> <X) limit of QCD is often one 
in which hard problems become tractable. However, the ground state of A^c = oo 
QCD is a chiral density wave, not a color superconductor |l9]. At asymptotically 
high densities color superconductivity persists up to Nc's of order thousands pd,bl| 
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nitude of the gaps in quark matter at large but accessible /i has been esti- 
mated using two broad strategies. The first class of estimates are done within 
the context of models whose parameters are chosen to reproduce zero density 
physics ^,^ f f |,^,^ 13 I^, p3| , p4|j2^ , p^ . The second strategy for estimating gaps 



and critical temperatures is to use /i = oo physics as a guide. At asymptot- 
ically large fj,, models with short-range interactions are bound to fail because 
the dominant interaction is due to the long-range magnetic interaction coming 
from single-gluon exchange |12| , |26| . The coUinear infrared divergence in small an- 
gle scattering via one-gluon exchange (which is regulated by dynamical screen- 
ing [g6|) results in a gap which is parametrically larger at /x — > oo than it 
would be for any point- like four-fcrmion interaction JEO] . Weak coupling esti- 



mates of the gap @,||,||J2§|3^,|l|,||j3|,||,|||6|3^,|§ are valid at asymptot 
ically high densities, with chemical potentials /i ^ 10^ MeV 37 1. Neither class 
of methods can be trusted quantitatively for quark number chemical potentials 
/x ~ 400 — 500 MeV, as appropriate for the quark matter which may occur in the 
cores of neutron stars. It is nevertheless satisfying that two very different ap- 
proaches, one using zero density phenomenology to normalize models, the other 
using weak-coupling methods valid at asymptotically high density, yield predic- 
tions for the gaps and critical temperatures at accessible densities which are in 
good agreement: the gaps at the Fermi surface are of order tens to 100 MeV, 
with critical temperatures about half as large. 

Tc ~ 50 MeV is much larger relative to the Fermi momentum than in low 
temperature superconductivity in metals. This reflects the fact that color super- 
conductivity is induced by an attraction due to the primary, strong, interaction 
in the theory, rather than having to rely on much weaker secondary interactions, 
as in phonon mediated superconductivity in metals. Quark matter is a high-Tc 
superconductor by any reasonable definition. Its Tc is nevertheless low enough 
that it is unlikely the phenomenon can be realized in heavy ion collisions. 



2 Color Superconductivity in Compact Stars 

Our current understanding of the color superconducting state of quark matter 
leads us to believe that it may occur naturally in compact stars. The critical tem- 
perature Tc below which quark matter is a color superconductor is high enough 
that any quark matter which occurs within neutron stars that are more than 
a few seconds old is in a color superconducting state. In the absence of lattice 
simulations, present theoretical methods are not accurate enough to determine 
whether neutron star cores are made of hadronic matter or quark matter. They 
also cannot determine whether any quark matter which arises will be in the CFL 
or 2SC phase: the difference between the u, d and s Fermi momenta will be a 



before being supplanted by the phase described in Ref. |ig|. At any finite Nc, color 
superconductivity occurs at arbitrarily weak coupling whereas the chiral density 
wave does not. For Nc = 3, color superconductivity is still favored over the chiral 
density wave (although not by much) even if the interaction is so strong that the 
color superconductivity gap is ^ fi/2 pd]. 
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few tens of MeV which is comparable to estimates of the gap A; the CFL phase 
occurs when A is large compared to all differences between Fermi momenta. Just 
as the higher temperature regions of the QCD phase diagram are being mapped 
out in heavy ion collisions, we need to learn how to use neutron star phenom- 
ena to determine whether they feature cores made of 2SC quark matter, CFL 
quark matter or hadronic matter, thus teaching us about the high density region 
of the QCD phase diagram. It is therefore important to look for astrophysical 
consequences of color superconductivity. 



2.1 Equation of State 

Much of the work on the consequences of quark matter within a compact star has 
focussed on the effects of quark matter on the equation of state, and hence on 
the radius of the star. As a Fermi surface phenomenon, color superconductivity 
has little effect on the equation of state: the pressure is an integral over the 
whole Fermi volume. Color superconductivity modifies the equation of state at 
the ~ {A/^Y level, typically by a few percent 1^. Such small effects can be 
neglected in present calculations, and for this reason we will not attempt to 
survey the many ways in which observations of neutron stars are being used to 
constrain the equation of state ||3^] . 

We will describe one current idea, however. As a neutron star in a low mass 
X-ray binary (LMXB) is spun up by accretion from its companion, it becomes 
more oblate and its central density decreases. If it contains a quark matter core, 
the volume fraction occupied by this core decreases, the star expands, and its 
moment of inertia increases. This raises the possibility pQ | of a period during the 
spin-up history of an LMXB when the neutron star is gaining angular momen- 
tum via accretion, but is gaining sufficient moment of inertia that its angular 
frequency is hardly increasing. In their modelling of this effect, Glendenning 
and Weber M] discover that LMXB's should spend a significant fraction of 
their history with a frequency of around 200 Hz, while their quark cores are 
being spun out of existence, before eventually spinning up to higher frequencies. 
This may explain the observation that LMXB frequencies are clustered around 



250-350 Hz 41 , which is otherwise puzzling in that it is thought that LMXB's 
provide the link between canonical pulsars and millisecond pulsars, which have 
frequencies as large as 600 Hz [Q. It will be interesting to see how robust the 
result of Ref. |Q is to changes in model assumptions and also how its predictions 
fare when compared to those of other explanations which posit upper bounds on 
LMXB frequencies |43| , rather than a most probable frequency range with no as- 
sociated upper bound |Q. We note here that because Glendenning and Weber's 
effect depends only on the equation of state and not on other properties of quark 
matter, the fact that the quark matter must in fact be a color superconductor 
will not affect the results in any significant way. If Glendenning and Weber's 
explanation for the observed clustering of LMXB frequencies proves robust, it 
would imply that pulsars with lower rotational frequencies feature quark matter 
cores. 
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2.2 Cooling by Neutrino Emission 

We turn now to neutron star phenomena which are affected by Fermi surface 
physics. For the first 10^^^ years of its life, the cooUng of a neutron star is 
governed by the balance between heat capacity and the loss of heat by neutrino 
emission. How are these quantities affected by the presence of a quark matter 
core? This has been addressed recently in Refs. Q,Q, following earlier work 
in Ref. |46|. Both the specific heat Cy and the neutrino emission rate L^ are 
dominated by physics within T of the Fermi surface. If, as in the CFL phase, 
all quarks have a gap Z\ ^ T then the contribution of quark quasiparticles to 
Cy and Li, is suppressed by ~ exp(— Z\/r). There may be other contributions 
to Li, |4J], but these are also very small. The specific heat is dominated by that 
of the electrons, although it may also receive a small contribution from the CFL 
phase Goldstone bosons. Although further work is required, it is already clear 
that both Cy and Li, are much smaller than in the nuclear matter outside the 
quark matter core. This means that the total heat capacity and the total neutrino 
emission rate (and hence the cooling rate) of a neutron star with a CFL core 
will be determined completely by the nuclear matter outside the core. The quark 
matter core is "inert" : with its small heat capacity and emission rate it has little 
influence on the temperature of the star as a whole. As the rest of the star emits 
neutrinos and cools, the core cools by conduction, because the electrons keep it 
in good thermal contact with the rest of the star. These qualitative expectations 
are nicely borne out in the calculations presented by Page et al. [ [f5[ . 

The analysis of the cooling history of a neutron star with a quark matter 
core in the 2SC phase is more complicated. The red and green up and down 
quarks pair with a gap many orders of magnitude larger than the temperature, 
which is of order 10 keV, and are therefore inert as described above. Any strange 
quarks present will form a (ss) condensate with angular momentum J = 1 which 
locks to color in such a way that rotational invariance is not broken p7|. The 



resulting gap has been estimated to be of order hundreds of keV 47 1, although 
applying results of Ref. Es] suggests a somewhat smaller gap, around 10 keV. 
The blue up and down quarks also pair, forming a J = 1 condensate which breaks 
rotational invariance S . The related gap was estimated to be a few keV H , but 
this estimate was not robust and should be revisited in light of more recent 
developments given its importance in the following. The critical temperature 
Tc above which no condensate forms is of order the zero-temperature gap A. 



[Tc — 0.57Z\ for J = condensates |27 .) Therefore, if there are quarks for 



which Z\ ^ T or smaller, these quarks do not pair at temperature T. Such 
quark quasiparticles will radiate neutrinos rapidly (via direct URCA reactions 
like d— >M + e + i?, M— >(i + e+ + i^, etc.) and the quark matter core will cool 
rapidly and determine the cooling history of the star as a whole [^^Q . The star 
will cool rapidly until its interior temperature is T < Tc ~ Z\, at which time 
the quark matter core will become inert and the further cooling history will be 
dominated by neutrino emission from the nuclear matter fraction of the star. If 
future data were to show that neutron stars first cool rapidly (direct URCA) and 
then cool more slowly, such data would allow an estimate of the smallest quark 
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matter gap. We are unlikely to be so lucky. The simple observation of rapid 
cooling would not be an unambiguous discovery of quark matter with small 
gaps; there are other circumstances in which the direct URCA processes occur. 
However, if as data on neutron star temperatures improves in coming years the 
standard cooling scenario proves correct, indicating the absence of the direct 
URCA processes, this would rule out the presence of quark matter with gaps in 
the 10 keV range or smaller. The presence of a quark matter core in which all 
gaps are ^ T can never be revealed by an analysis of the cooling history. 

2.3 Supernova Neutrinos 

We now turn from neutrino emission from a neutron star which is many years 
old to that from the protoncutron star during the first seconds of a supernova. 
Carter and Reddy Efl] have pointed out that when this protoncutron star is 
at its maximum temperature of order 30-50 MeV, it may have a quark mat- 
ter core which is too hot for color superconductivity. As such a protoneutron 
star core cools over the next few seconds, this quark matter will cool through 
Tc, entering the color superconducting regime of the QCD phase diagram. For 
T r^ Tc, the specific heat rises and the cooling slows. Then, as T drops fur- 
ther and A increases to become greater than T, the specific heat drops rapidly. 
Furthermore, as the number density of quark quasiparticles becomes suppressed 
by exp(— Z\/T), the neutrino transport mean free path rapidly becomes very 
long 1)4 9| . This means that all the neutrinos previously trapped in the now color 
superconducting core are able to escape in a sudden burst. If a terrestrial neu- 
trino detector sees thousands of neutrinos from a future supernova. Carter and 
Reddy's results suggest that there may be a signature of the transition to color 
superconductivity present in the time distribution of these neutrinos. Neutrinos 
from the core of the protoneutron star will lose energy as they scatter on their 
way out, but because they will be the last to reach the surface of last scattering, 
they will be the final neutrinos received at the earth. If they are released from 
the quark matter core in a sudden burst, they may therefore result in a bump at 
late times in the temporal distribution of the detected neutrinos. More detailed 
study remains to be done in order to understand how Carter and Reddy's sig- 
nature, dramatic when the neutrinos escape from the core, is processed as the 
neutrinos traverse the rest of the protoncutron star and reach their surface of 
last scattering. 

2.4 R-mode Instabilities 

Another arena in which color superconductivity comes into play is the physics of 
r-mode instabilities. A neutron star whose angular rotation frequency fi is large 
enough is unstable to the growth of r-mode oscillations which radiate away angu- 
lar momentum via gravitational waves, reducing f2. What does "large enough" 
mean? The answer depends on the damping mechanisms which act to prevent 
the growth of the relevant modes. Both shear viscosity and bulk viscosity act to 
damp the r-modes, preventing them from going unstable. The bulk viscosity and 
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the quark contribution to the shear viscosity both become exponentially small in 
quark matter with A > T and as a result, as Madsen pO| has shown, a compact 
star made entirely of quark matter with gaps A — 1 MeV or greater is unstable if 
its spin frequency is greater than tens to 100 Hz. Many compact stars spin faster 
than this, and Madsen therefore argues that compact stars cannot be strange 
quark stars unless some quarks remain ungapped. Alas, this powerful argument 
becomes much less powerful in the context of a neutron star with a quark matter 
core. First, the r-mode oscillations have a wave form whose amplitude is largest 
at large radius, outside the core. Second, in an ordinary neutron star there is a 
new source of damping: friction at the boundary between the crust and the neu- 
tron superfluid "mantle" keeps the r-modes stable regardless of the properties 
of a quark matter core |5^,Q . 

2.5 Magnetic Field Evolution 

Next, we turn to the physics of magnetic fields within color superconducting neu- 
tron star cores [|2|J5^ . The interior of a conventional neutron star is a superfluid 
(because of neutron-neutron pairing) and is an electromagnetic superconductor 
(because of proton- proton pairing). Ordinary magnetic fields penetrate it only 
in the cores of magnetic flux tubes. A color superconductor behaves diff'erently. 
At first glance, it seems that because a diquark Cooper pair has nonzero elec- 
tric charge, a diquark condensate must exhibit the standard Meissner effect, 
expelling ordinary magnetic fields or restricting them to flux tubes within whose 
cores the condensate vanishes. This is not the case p3|. In both the 2SC and 
CFL phase a linear combination of the f/(l) gauge transformation of ordinary 
electromagnetism and one (the eighth) color gauge transformation remains un- 
broken even in the presence of the condensate. This means that the ordinary 
photon Afj^ and the eighth gluon G^ are replaced by new linear combinations 

A^ = cosao Af^ + sinao G^ 

A^ =- sin ao A^ + cos ao Gl (2) 

where A^ is massless and Arf is massive. That is, Bq satisfies the ordinary 
Maxwell equations while Bx experiences a Meissner effect. The mixing angle ao 
is the analogue of the Weinberg angle in electroweak theory, in which the pres- 
ence of the Higgs condensate causes the A^ and the third SU{2)iy gauge boson 
to mix to form the photon, A^, and the massive Z boson, sin(ao) is proportional 
to e/g and turns out to be about 1/20 in the 2SC phase and 1/40 in the CFL 
phase p3[ . This means that the Q'Photon which propagates in color supercon- 
ducting quark matter is mostly photon with only a small gluon admixture. If 
a color superconducting neutron star core is subjected to an ordinary magnetic 
field, it will either expel the X component of the flux or restrict it to flux tubes, 
but it can (and does |£3|) admit the great majority of the flux in the form of 
a Bq magnetic field satisfying Maxwell's equations. The decay in time of this 
"free field" (i.e. not in flux tubes) is limited by the Q-conductivity of the quark 
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matter. A color superconductor is not a Q-superconductor — that is the whole 
point — but it turns out to be a very good Q-conductor due to the presence 
of electrons: the Bq magnetic field decays only on a time scale which is much 
longer than the age of the universe |53] . This means that a quark matter core 
within a neutron star serves as an "anchor" for the magnetic field: whereas in 
ordinary nuclear matter the magnetic flux tubes can be dragged outward by the 
neutron superfluid vortices as the star spins down p3 , the magnetic flux within 
the color superconducting core simply cannot decay. Even though this distinc- 
tion is a qualitative one, it will be difficult to confront it with data since what is 
observed is the total dipole moment of the neutron star. A color superconducting 
core anchors those magnetic flux lines which pass through the core, while in a 
neutron star with no quark matter core the entire internal magnetic field can 
decay over time. In both cases, however, the total dipole moment can change 
since the magnetic flux lines which do not pass through the core can move. 



3 Crystalline Color Superconductivity and Glitches in 
Quark Matter 

The final consequence of color superconductivity we wish to discuss is the possi- 
bility that (some) glitches may originate within quark matter regions of a com- 
pact star pq] . In any context in which color superconductivity arises in nature, it 
is likely to involve pairing between species of quarks with differing chemical po- 
tentials. If the chemical potential difference is small enough, BCS pairing occurs 
as we have been discussing. If the Fermi surfaces are too far apart, no pairing 
between the species is possible. The transition between the BCS and unpaired 
states as the splitting between Fermi momenta increases has been studied in 
electron superconductors [ p5[ , nuclear superfluids [ p6[ and QCD superconduc- 
tors [plp,M, assuming that no other state intervenes. However, there is good 
reason to think that another state can occur. This is the "LOFF" state, first 
explored by Larkin and Ovchinnikov [b8| and Fulde and Ferrell pS] in the con- 
text of electron superconductivity in the presence of magnetic impurities. They 
found that near the unpairing transition, it is favorable to form a state in which 
the Cooper pairs have nonzero momentum. This is favored because it gives rise 
to a region of phase space where each of the two quarks in a pair can be close 
to its Fermi surface, and such pairs can be created at low cost in free energy. 
Condensates of this sort spontaneously break translational and rotational invari- 
ance, leading to gaps which vary periodically in a crystalline pattern. If in some 
shell within the quark matter core of a neutron star (or within a strange quark 
star) the quark number densities are such that crystalline color superconductiv- 
ity arises, rotational vortices may be pinned in this shell, making it a locus for 
glitch phenomena. 

In Ref . iQ , we have explored the range of parameters for which crystalline 
color superconductivity occurs in the QCD phase diagram, upon making vari- 
ous simplifying assumptions. We focus primarily on a toy model in which the 
quarks interact via a four-fcrmion interaction with the quantum numbers of sin- 
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gle gluon exchange. Also, we only consider pairing between u and d quarks, with 
l^d = P- + ^IJ- and /i„ = yU — Sjj,, whereas we expect a LOFF state wherever the 
difference between the Fermi momenta of any two quark flavors is near an un- 
pairing transition, including, for example, near the unlocking phase transition 
between the 2SC and CFL phases. 

In the LOFF state, each Cooper pair carries momentum 2q with |q| « 1.25/x. 
The condensate and gap parameter vary in space with wavelength 7r/|q|. In 
Ref. Q , we simplify the calculation by assuming that the condensate varies in 
space like a plane wave, leaving the determination of the crystal structure of 
the QCD LOFF phase to future work. We give an ansatz for the LOFF wave 
function, and by variation obtain a gap equation which allows us to solve for the 
gap parameter Aa, the free energy and the values of the diquark condensates 
which characterize the LOFF state at a given 6fj, and |q|. We then vary |q|, to 
find the preferred (lowest free energy) LOFF state at a given Sfj,, and compare the 
free energy of the LOFF state to that of the BCS state with which it competes. 
We show results for one choice of parameters^ in Fig. 1(a). The LOFF state 
is characterized by a gap parameter Aa and a diquark condensate, but not by 
an energy gap in the dispersion relation: we obtain the quasiparticle dispersion 
relations p8[ and find that they vary with the direction of the momentum, 
yielding gaps that vary from zero up to a maximum of A a- The condensate is 
dominated by the regions in momentum space in which a quark pair with total 
momentum 2q has both members of the pair within '^ Aa of their respective 
Fermi surfaces. 

Because it violates rotational invariance by involving Cooper pairs whose 
momenta are not antiparallel, the quark matter LOFF state necessarily features 
condensates in both the J — and J — I channels. (Cooper pairs in the sym- 
metric J = 1 channel are antisymmetric in color but symmetric in flavor, and 
are impossible in the original LOFF context of pairing between electrons, which 
have neither color nor flavor.) Both J — and J = 1 condensates are present 
even if there is no interaction in the J = 1 channel, as is the case when we 
use a four-fermion interaction with the quantum numbers of Lorentz-invariant 
single gluon exchange. Because there is no interaction in the J — 1 channel, the 
J = 1 condensate does not affect the quasiparticle dispersion relations; that is, 
the J = 1 gap parameter vanishes. 

The LOFF state is favored for values of (5/i which satisfy (5/ii < Sfj. < (5/i2 
as shown in Fig. 1(b), with 5/zi/Z\o = 0.707 and 6112/AQ — 0.754 in the weak 
coupling limit in which Aq <C fJ-- (For Sfi < Sfii, we have the 2SC phase with 
gap Aq.) At weak coupling, the LOFF gap parameter decreases from 0.23Z\o at 
Sfi ~ (5/ii (where there is a first order BCS-LOFF phase transition) to zero at 



^ Our model Hamiltonian has two parameters, the four-fermion coupling G and a 
cutoff A. We often use the value of Aq, the BCS gap obtained at Sfi = 0, to describe 
the strength of the interaction: small Ao corresponds to small G. When we wish to 
study the dependence on the cutoff, we vary A while at the same time varying the 
coupling G such that Ao is kept fixed. We expect that the relation between other 
physical quantities and Aq will be reasonably insensitive to variation of A. 
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Fig. 1. (a) LOFF and BCS gap parameters as a function of S^, with coupling chosen 
so that Ao = 40 MeV. At each 5fi we have varied |q| to find the LOFF state with the 
lowest free energy. The vertical dashed line marks Sn = Sfii, the value of Sfi above 
which the LOFF state has lower free energy than BCS. (b) The interval of Sfi within 
which the LOFF state occurs as a function of the coupling, parametrized by the BCS 
gap Ao shown in GeV. Here and in (a), the average quark chemical potential p, has 
been set to 0.4 GeV, correspondin g to a baryon density of about 4 to 5 times that in 
nuclear matter. A crude estimate |48| suggests that in quark matter at this density, 
5/x ~ 15 — 30 MeV depending on the value of the density-dependent effective strange 
quark mass. Below the solid line, there is a LOFF state. Below the dashed line, the 
BCS state is favored. The different lines of each type correspond to different cutoffs on 
the four-fermion interaction: A = 0.8 GeV to 1.6 GeV. 5fii/Ao and 8^2/ Ao show little 
cutoff-dependence, and the cutoff-dependence disappears completely as Ao,5yL -^ 0. 
A = 1 GeV in (a). 



(5/i = 5/i2 (where there is a second order LOFF- normal transition). Except for 
very close to (5/i2 , the critical temperature above which the LOFF state melts will 
be much higher than typical neutron star temperatures. At stronger coupling the 
LOFF gap parameter decreases relative to Z\o and the window of Sji/Aq within 
which the LOFF state is favored shrinks, as seen in Fig. 1(b). The single gluon 
exchange interaction used in Fig. 1 is neither attractive nor repulsive in the J = 1 
channel: the width of the LOFF window grows if the interaction is modified to 
include an attraction in this channel [Esl . 

Near the second-order critical point 5^2 , we can describe the phase transition 
with a Ginzburg-Landau effective potential. The order parameter for the LOFF- 
to-normal phase transition is 



1>{r) = -i(e.6e„/33^'^"(r)C75^'"'(r)) 



(3) 



so that in the normal phase ^(r) = 0, while in the LOFF phase <?(r) = Fag^^'^'^ . 
(The gap parameter is related to the order parameter by A a = GFa-) Expressing 
the order parameter in terms of its Fourier modes <^(k), we write the LOFF free 
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energy (relative to the normal state) as 

Fi{3{k)}) = j2 {c2ik')\m\' + c.iemkr + oim) . (4) 

k 

For (5/i > (5/i2, all of the C2(fc^) are positive and the normal state is stable. Just 
below the critical point, all of the modes ^(k) are stable except those on the 
sphere |k| = 2q2, where q2 is the value of |q| at Sfj.2 (so that 52 — l-25/i2 — 0.9Z\o 
at weak coupling). In general, many modes on this sphere can become nonzero, 
giving a condensate with a complex crystal structure. We consider the simplest 
case of a plane wave condensate where only the one mode ^(k — 2q2) = Fa is 
nonvanishing. Dropping all other modes, we have 

F{rA) = a{6fi - 6fi2){FAf + b{FAf (5) 

where a and h are positive constants. Finding the minimum-energy solution for 
djjL < (5/X2, we obtain simple power-law relations for the condensate and the free 
energy: 

FA{Sti)^Kr{St^2-6fi)^/^, F{d^i) = -KF{Sfi2-Sti)^. (6) 

These expressions agree well with the numerical results obtained by solving the 
gap equation Q. The Ginzburg-Landau method does not specify the propor- 
tionality factors Kr and Kp, but analytical expressions for these coefficients 
can be obtained in the weak coupling limit by explicitly solving the gap equa- 
tion MM, yielding 



Kf = {4fiy7T^){{q2/5fi2? - 1) ^ 0.178/22. ^ ' 

Notice that because {Sfi2 — Sjii)/ 5^2 is small, the power-law relations (ph are a 
good model of the system throughout the entire LOFF interval 5^i < (5/i < 8^2 
where the LOFF phase is favored over the BCS phase. The Ginzburg-Landau 
expression (m gives the free energy of the LOFF phase near (5/i2 , but it cannot be 
used to determine the location 5^i of the first-order phase transition where the 
LOFF window terminates. (Locating the first-order point requires a comparison 
of LOFF and BCS free energies.) 

The quark matter which may be present within a compact star will be in 
the crystalline color superconductor (LOFF) state if (5/i/Z\o is in the requisite 
range. For a reasonable value of (5/i, say 25 MeV, this occurs if the gap Z\o which 
characterizes the uniform color superconductor present at smaller values of Sfj, 
is about 40 MeV. This is in the middle of the range of present estimates. Both 
6^ and Aq vary as a function of density and hence as a function of radius in 
a compact star. Although it is too early to make quantitative predictions, the 
numbers are such that crystalline color superconducting quark matter may very 
well occur in a range of radii within a compact star. It is therefore worthwhile 
to consider the consequences. 
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Many pulsars have been observed to glitch. Glitches are sudden jumps in 
rotation frequency il which may be as large as Afi/fi ~ 10~^, but may also 
be several orders of magnitude smaller. The frequency of observed glitches is 
statistically consistent with the hypothesis that all radio pulsars experience 



glitches 1 61 . Glitches are thought to originate from interactions between the 
rigid neutron star crust, typically somewhat more than a kilometer thick, and 
rotational vortices in a neutron superfluid. The inner kilometer of crust con- 
sists of a crystal lattice of nuclei immersed in a neutron superfluid |62|] . Because 
the pulsar is spinning, the neutron superfluid (both within the inner crust and 
deeper inside the star) is threaded with a regular array of rotational vortices. As 
the pulsar's spin gradually slows, these vortices must gradually move outwards 
since the rotation frequency of a superfluid is proportional to the density of vor- 
tices. Deep within the star, the vortices are free to move outwards. In the crust, 
however, the vortices are pinned by their interaction with the nuclear lattice. 
Models Q differ in important respects as to how the stress associated with 
pinned vortices is released in a glitch: for example, the vortices may break and 
rearrange the crust, or a cluster of vortices may suddenly overcome the pinning 
force and move macroscopically outward, with the sudden decrease in the angu- 
lar momentum of the superfluid within the crust resulting in a sudden increase 
in angular momentum of the rigid crust itself and hence a glitch. All the models 
agree that the fundamental requirements are the presence of rotational vortices 
in a superfluid and the presence of a rigid structure which impedes the motion of 
vortices and which encompasses enough of the volume of the pulsar to contribute 
significantly to the total moment of inertia. 



Although it is premature to draw quantitative conclusions, it is interesting to 
speculate that some glitches may originate deep within a pulsar which features 
a quark matter core, in a region of that core which is in a LOFF crystalline 
color superconductor phase. A three flavor analysis is required to estimate over 
what range of densities LOFF phases may arise, as either (ud), (us) or (ds) 
condensates approach their unpairing transitions. Comparison to existing mod- 
els which describe how p'^, pp and pp vary within a quark matter core in a 
neutron star [ p4[ would then permit an estimate of how much the LOFF region 
contributes to the moment of inertia of the pulsar. Furthermore, a three flavor 
analysis is required to determine whether the LOFF phase is a superfluid. If the 
only pairing is between u and d quarks, this 2SC phase is not a superfluid 0j8|, 
whereas if all three quarks pair in some way, a superfluid is obtained ||5|,|8[ . Hence- 
forth, we suppose that the LOFF phase is a superfluid, which means that if it 
occurs within a pulsar it will be threaded by an array of rotational vortices. It 
is reasonable to expect that these vortices will be pinned in a LOFF crystal, 
in which the diquark condensate varies periodically in space. Indeed, one of the 
suggestions for how to look for a LOFF phase in terrestrial electron supercon- 
ductors relies on the fact that the pinning of magnetic flux tubes (which, like the 
rotational vortices of interest to us, have normal cores) is expected to be much 
stronger in a LOFF phase than in a uniform BCS superconductor |6a|. 
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A real calculation of the pinning force experienced by a vortex in a crystalline 
color superconductor must await the determination of the crystal structure of 
the LOFF phase. We can, however, attempt an order of magnitude estimate 
along the same lines as that done by Anderson and Itoh IM for neutron vortices 
in the inner crust of a neutron star. In that context, this estimate has since 
been made quantitative | |67| , p8p^ . For one specific choice of parameters |48| , the 
LOFF phase is favored over the normal phase by a free energy Floff ^ 5 x 
(10 MeV)'* and the spacing between nodes in the LOFF crystal isb = 7r/(2|q|) ~ 
9 fm. The thickness of a rotational vortex is given by the correlation length 
^ ^ 1/A ^ 25 fm. The pinning energy is the difference between the energy of a 
section of vortex of length b which is centered on a node of the LOFF crystal 
vs. one which is centered on a maximum of the LOFF crystal. It is of order 
Ep ^ -Floff fo"^ ^ 4 MeV. The resulting pinning force per unit length of vortex 
is of order fp ~ Ep/b'^ ^ (4 MeV)/(80 fm^). A complete calculation will be 
challenging because b < £^, and is likely to yield an fp which is somewhat less 
than that we have obtained by dimensional analysis. Note that our estimate 
of fp is quite uncertain both because it is only based on dimensional analysis 
and because the values of Z\, 6 and Floff are uncertain. (We have a good 
understanding of all the ratios A/Aq, 6fi/Ao, q/Ao and consequently bAo in the 
LOFF phase. It is of course the value of the BCS gap Aq which is uncertain.) It 
is premature to compare our crude result to the results of serious calculations of 
the pinning of crustal neutron vortices as in Refs. p7| , p8|p^ . It is nevertheless 
remarkable that they prove to be similar: the pinning energy of neutron vortices 
in the inner crust is Fp « 1 — 3 MeV and the pinning force per unit length is 
/p « (1 - 3 MeV)/(200 - 400 fm^). 

The reader may be concerned that a glitch deep within the quark matter core 
of a neutron star may not be observable: the vortices within the crystalline color 
superconductor region suddenly unpin and leap outward; this loss of angular 
momentum is compensated by a gain in angular momentum of the layer outside 
the LOFF region; how quickly, then, does this increase in angular momentum 
manifest itself at the surface of the star as a glitch? The important point here 
is that the rotation of any supcrfluid region within which the vortices are able 
to move freely is coupled to the rotation of the outer crust on very short time 
scales [p3. This rapid coupling, due to electron scattering off vortices and the 
fact that the electron fluid penetrates throughout the star, is usually invoked to 
explain that the core nucleon supcrfluid speeds up quickly after a crustal glitch: 
the only long relaxation time is that of the vortices within the inner crust | |69| . 
Here, we invoke it to explain that the outer crust speeds up rapidly after a LOFF 
glitch has accelerated the quark matter at the base of the nucleon supcrfluid. 
After a glitch in the LOFF region, the only long relaxation times are those of 
the vortices in the LOFF region and in the inner crust. 

A quantitative theory of glitches originating within quark matter in a LOFF 
phase must await further calculations, in particular a three flavor analysis and 
the determination of the crystal structure of the QCD LOFF phase. However, 
our rough estimate of the pinning force on rotational vortices in a LOFF region 
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suggests that this force may be comparable to that on vortices in the inner crust 
of a conventional neutron star. Perhaps, therefore, glitches occurring in a region 
of crystalline color superconducting quark matter may yield similar phenomenol- 
ogy to those occurring in the inner crust. This is surely strong motivation for 
further investigation. 

Perhaps the most interesting consequence of these speculations arises in the 
context of compact stars made entirely of strange quark matter. The work of 
Witten |70J and Farhi and Jaffe |7^ raised the possibility that strange quark 
matter may be stable relative to nuclear matter even at zero pressure. If this 
is the case it raises the question whether observed compact stars — pulsars, for 
example — are strange quark stars [[72y73tl rather than neutron stars. A conven- 
tional neutron star may feature a core made of strange quark matter, as we have 
been discussing above.^ Strange quark stars, on the other hand, are made (al- 
most) entirely of quark matter with either no hadronic matter content at all or 
with a thin crust, of order one hundred meters thick, which contains no neutron 
superfluid |73| , |74| . The nuclei in this thin crust are supported above the quark 
matter by electrostatic forces; these forces cannot support a neutron fluid. Be- 
cause of the absence of superfluid neutrons, and because of the thinness of the 
crust, no successful models of glitches in the crust of a strange quark star have 
been proposed. Since pulsars are observed to glitch, the apparent lack of a glitch 
mechanism for strange quark stars has been the strongest argument that pulsars 
cannot be strange quark stars ||75|,|76 77 . This conclusion must now be revisited. 



Madsen's conclusion |Q that a strange quark star is prone to r-mode insta- 
bility due to the absence of damping must also be revisited, since the relevant 
oscillations may be damped within or at the boundary of a crystalline color 
superconductor region. 

The quark matter in a strange quark star, should one exist, would be a 
color superconductor. Depending on the mass of the star, the quark number 
densities increase by a factor of about two to ten in going from the surface to the 
center ]73] . This means that the chemical potential differences among the three 
quarks will vary also, and there could be a range of radii within which the quark 
matter is in a crystalline color superconductor phase. This raises the possibility 
of glitches in strange quark stars. Because the variation in density with radius 
is gradual, if a shell of LOFF quark matter exists it need not be particularly 
thin. And, we have seen, the pinning forces may be comparable in magnitude 
to those in the inner crust of a conventional neutron star. It has recently been 
suggested (for reasons unrelated to our considerations) that certain accreting 
compact stars may be strange quark stars |7^ , although the evidence is far from 
unambiguous ||7g|. In contrast, it has been thought that, because they glitch, 
conventional radio pulsars cannot be strange quark stars. Our work questions 



Note that a convincing discovery of a quark matter core within an otherwise hadronic 
neutron star would demonstrate conclusively that strange quark matter is not stable 
at zero pressure, thus ruling out the existence of strange quark stars. It is not possible 
for neutron stars with quark matter cores and strange quark stars to both be stable. 



16 



Mark Alford, Jeffrey A. Bowers, and Krishna Rajagopal 



this assertion by raising the possibihty that ghtches may originate within a layer 
of quark matter which is in a crystaUine color superconducting state. 

There has been much recent progress in our understanding of how the pres- 
ence of color superconducting quark matter in a compact star would affect five 
different phenomena: cooling by neutrino emission, the pattern of the arrival 
times of supernova neutrinos, the evolution of neutron star magnetic fields, r- 
mode instabilities and glitches. Nevertheless, much theoretical work remains to 
be done before we can make sharp proposals for which astrophysical observations 
can teach us whether compact stars contain quark matter, and if so whether it 
is in the 2SC or CFL phase. 
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